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A  PROPOSED  STATIONARY  OR  MOBILE 
OIL  SHALE  RETORTING  SYSTEM 


Abstract 


An  investigation  of  the  minimum  power  required  for  retorting  oil  shale  in  situ 
using  the  thermal  heat  from  a  nuclear  reactor  revealed  the  necessity  of  having  a  340  MW 
reactor.  In  order  to  map  temperature  profiles  in  a  nuclear  chimney,  a  new  equation 
was  derived  which  computes  the  ratio  of  the  solid -to -fluid  temperature.  Heat  exchangers 
and  pumping  power  requirements  were  established  for  the  model  under  consideration. 

The  cost  of  the  nuclear  retorting  system  is  negligible  when  prorated  over  the  lifetime 
of  the  reactor. 


Introduction 


Earlier  studies  at  the  Lawrence  Radiation  Laboratory  have  shown  that  it  is 
feasible  to  extract  oil  from  the  broken  oil  shale  rock  in  a  chimney  created  by  an  under¬ 
ground  nuclear  explosion.^  One  method  suggests  that  a  significant  amount  of  oil  from 

large  chunks  of  oil  shale  can  be  removed  by  igniting  the  broken  shale  and  maintaining 
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combustion  with  air  supplied  to  the  chimney  through  drill  holes.  By  this  process, 
called  retorting,  the  heat  from  the  combustion  process  would  cause  the  solid  hydro¬ 
carbon  in  the  shale  to  melt,  producing  oil  which  flows  by  gravity  to  the  base  of  the 
chimney  where  it  will  be  drained  and  pumped  to  surface  storage  facilities.  Other 
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techniques  for  retorting  oil  shale  have  been  reported  where  hot  gases,  steam,  and 
electrical  heaters  were  employed. 

Instead  of  utilizing  heat  from  a  combustion-driven  front  to  enhance  the  retorting 
process  or  steam  which  may  cause  undue  chemical  reactions,  oil  from  the  field  may  be 
utilized  in  a  heat  exchange  process  with  a  nuclear  reactor  to  liberate  encapsulated  oil 
in  oil  shale.  An  arbitrary  retorting  temperature  of  700°F  was  chosen  after  an  examina- 
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tion  of  experimental  data  reviewed  by  Williamson.  Some  investigators  found  that  no 
oil  was  formed  from  kerogen  when  shales  were  heated  while  mixed  with  carnotite: 
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(1)  for  28  days  at  392  to  644®F 

(2)  for  26  days  at  176°F 

(3)  for  30  hr  at  500®F  with  various  catalysts 

(4)  for  one  week  at  509®F. 

The  following  description  of  the  general  reaction  of  kerogen  to  temperature 
increases  was  reported: 

617°F  Kerogen  softens,  swells,  becomes  black,  has  general 

appearance  of  soft  bitumen.  Solid  nitrogen  compounds  are 
converted  to  soluble  compounds. 

572  —  761°F  The  "water  producing"  period — water-of-combination  and 
ammonia  are  released.  Easily  volatilized  organic  nitrogen 
compounds  and  hydrogen  sulfide  are  formed. 

734°F  Heavy  oil  vapors  and  gas  begin  to  be  formed  and  distilled. 

932  —  1112°F  Production  of  condensible  oil  vapors  ceases. 

+  1112®F  Noncondensible  hydrocarbons  and  ammonia  continue  to  be 

produced. 

It  was  also  reported  that  for  the  rapid  retorting  of  Colorado  shales,  oil  started 
to  be  expelled  from  the  shale  at  689°F  and  ceased  at  977°F.  For  slow  retorting  of 
shale  from  the  Colorado  fields,  oil  began  to  be  expelled  at  680®F  and  ceased  at  860°F. 

No  specific  times  were  reported  which  distinguished  slow  or  fast  retorting  but  retorting 
times  vary  as  the  square  of  particle  diameter. 

This  report  proposes  a  system  whereby  readily  available  oil  may  be  passed 
through  a  heat  exchanger  coupled  to  a  nuclear  reactor,  cycled  through  a  nuclear 
explosion-derived  oil  shale  chimney  to  liberate  oil  from  the  shale. 

The  merit  of  this  system  is  in  the  availability  of  most  of  its  component  parts  as 
off-the-shelf  items.  None  of  the  items  indicated  in  the  schematic  (Fig.  1)  require  an 
extensive  research  and  development  phase;  only  an  engineering  design  phase  is  required. 
The  nuclear  reactor  and  heat  exchangers  are  items  similar  to  those  developed  by  the 
AEC  where  liquid  sodium  was  the  primary  reactor  heat  exchange  medium  and  graphite 
the  moderator  material. 

There  are  three  heat  exchange  passes  in  this  system.  The  primary  cycle  of 
liquid  sodium  which  removes  heat  from  the  reactor  may  become  radioactive,  the 
secondary  cycle  of  sodium-potassium  (Na-K)  eliminates  the  possibility  of  radioactive 
contamination  and  hazard,  and  the  third  cycle  or  oil  heating  cycle  transfers  the  reactor 
energy  to  the  oil  shale. 

The  nuclear  reactor  and  heat  exchanger  units  should  be  as  compact  as  possible 
and  fully  shielded.  A  single  reactor  retorting  unit  may  be  sufficient  for  the  total  oil 
shale  field  or  for  isolated  oil  shale  formations  that  show  an  economic  promise. 
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Nuclear  reactor 


Heat  exchanger 


Method  of  Operation 


The  retorting  medium  would  be  oil 
(hydrocarbon)  readily  derived  from  on- 
the-site  natural  oil  operations  or  externally 
supplied.  Residual  heat  from  the  effects 
of  a  nuclear  device  may  cause  partial 
retorting.  Oil  accumulation  from  adjacent 
wells  may  be  sufficient  to  start  the  con¬ 
trolled  retorting  process.  Otherwise, 
oil  would  be  stored  in  a  reservoir  of 


^  ^  Heat  exchanger  i  sufficient  quantity  (as  shown  in  Fig.  2) 

I  to  perform  the  required  heat  transporting 

.  gipver^low  function  of  raising  the  temperature  of 

reservoir  6 

L____J  the  fragmented  shale  to  cause  the  con- 

tained  kerogen  to  become  fluid.  The 

'  A  heating  of  the  oil  would  be  performed  by 

the  heat  exchanger  which  derives  its 
energy  from  the  secondary  loop  of  the 

Nuclear 

chimney  nuclear  circuit.  A  pump  would  draw  the 

oil  from  the  oil  storage  reservoir 
through  the  heat  exchanger  where  its 

w  temperature  would  be  raised  to  a 

desirable  level  and  would  inject  the  hot 

_ ^ _  oil  into  the  shale  pile  in  a  random  or 

^  Pump  spray  manner.  The  hot  oil  would  flow 

through  the  void  spaces  of  the  nuclear 
shattered  oil  shale,  raising  the  temperature 
Fig.  1.  Mobile  oil  shale  retorting  system.  of  the  shale  deposits  as  it  flows.  The  oil 

initially  injected  would  suffer  a  tempera¬ 
ture  drop  in  giving  up  its  heat  to  the  oil  shale  and  the  surrounding  walls.  A  collector 
pump  would  be  located  such  that  oil  which  had  passed  through  the  flooded  shale  bed 
would  be  returned  to  the  oil  reservoir.  The  oil  reheat  cycle  would  begin  again.  The 
continuous  injection  of  hot  oil  would  soon  raise  the  temperature  of  the  total  oil  shale 
bed  to  the  point  where  sufficient  fluidity  would  be  exhibited  by  the  oil  in  the  shale  mass. 
Excess  oil  or  oil  derived  from  the  shale  mass  would  be  pumped  as  overflow  to  a 
processing  plant  or  new  storage  bins.  After  the  retorted  oil  has  been  depleted  from  the 
nuclear  chimney,  processing  may  begin  immediately  on  a  new  chimney.  For  optimum 
operation  the  retorting  system  should  be  strategically  placed  within  an  array  of  nuclear 
chimneys.  The  nuclear  reactor  would  operate  continuously  over  the  lifetime  of  the  system. 


reservoir 


Overflow 
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Nuclear 

chimney 


Fig.  1.  Mobile  oil  shale  retorting  system. 
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Inlet 


Fig.  2.  Oil  shale  retorting  by  means  of  a  nuclear  reactor. 
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Chimney  Model 


The  size  and  shape  of  nuclear-explosive-derived  cavities  are  inferred  from  data 
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reported  at  the  Third  Plowshare  Symposium.  The  equation  for  scaling  cavity  radii  is 
reported  to  be 


R 


C 


wl/3 

(ph)^ 


(1) 


where 

R  =  cavity  radius  (m) 

W  =  yield  (kt) 

p  =  average  of  overburden  density  (g/cm  ) 
h  =  depth  of  burial  of  nuclear  device  (m) 

C  =  a  constant  to  be  empirically  derived. 

The  chimney  appears  to  be  of  cylindrical  shape  with  a  radius  approximately  equal  to 
the  initial  nuclear-explosive-derived  cavity.  The  chimney  height  is  approximated  by 


H  =  KR  , 


where  K  is  empirically  derived.  If  the  volume  of  the  cylinder  is  taken  to  be 


V  =  ; 


then 


y  = 


KttC^W 


(ph) 


374" 


(2) 


If  <j)  is  the  percent  porosity  of  the  nuclear  chimney,  then  the  average  void  volume  is 


V^  =  <t>V. 

A  100-kt  nuclear  device  would  create  a  nuclear  chimney  with  the  approximate  dimensions 
of 

R  =  43  m  (140  ft) 

H  =  221  m  (725  ft). 
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This  represents  a  volume  of  4.4  X  10  ft  of  oil  shale.  If  the  chimney  has  an  average 

9 

porosity  of  27%,  then  the  approximate  oil  shale  weight  to  be  retorted  is  5  X  10  lb. 

In  the  calculations  being  presented,  it  is  assumed  that  the  oil  yield  of  the  shale 
is  28  (gpt). 


Heat  Transfer  Calculations 


The  parameters  in  Table  I  are  assumed  in  the  oil  and  oil  shale  heat  transfer 
calculations  that  follow.  A  list  of  symbols  for  the  entire  report  appears  in  Appendix  A. 


Table  I.  Oil  and  shale  thermodynamic  properties. 


Parameter 

Unit 

Colorado 
shale  oil6 

Mid-Continental 

oil"^ 

Colorado 
oil  shale8,9,10 

Boiling  point^ 

op 

734 

900 

— 

Characterization 

factor 

11.87 

12.2 

— 

Density 

Ib/ft^ 

55.687 

54,938 

144 

Enthalpy 

Btu/lb 

— 

— 

232 

Gravity 

°API 

27 

29 

— 

Group 

III^ 

II 

— 

Melting  point 

op 

680 

— 

— 

Molecular  weight 

304 

310 

— 

Specific  heat 

Btu/lb  “F 

0.765/0.675^ 

0.775/0.675^ 

0.313  (700°F) 

Temperature 

op 

7  00 

700 

80 

Thermal 

conductivity 

Btu/ft  hr  °F 

0.0794 

0.0633 

1.23 

Vapor  pressure 

atm 

0.7 

0.1 

— 

Viscosity 

Ib/hr  ft 

0.6776 

0.4114 

— 

^Data  are  based  on  gas  oils  and  heavy  crude  oils  which  indicate  boiling  points  in 
excess  of  500°F. 

^Has  the  same  characteristics  as  Group  III. 

^Liquid /vapor. 

The  objective  of  this  heat  transfer  analysis  of  hot  oil  passing  through  a  porous 
oil  shale  bed  is  to  establish  the  minimum  power  requirements  for  a  nuclear  reactor  to 
supply  heat  energy  for  retorting  the  oil  shale. 

By  using  techniques  developed  by  Zabrodsky^^  fluid  flow  and  heat  transfer 
characteristics  will  be  determined.  Although  Zabrodsky’s  approach  is  based  on  the 
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theory  of  fluidized  beds,  his  formulations  assume  a  static  balance  between  the  heating 

medium  and  the  bed  particles;  i.  e. ,  the  fluid  temperature  leaving  the  bed  on  its  initial 

pass  is  assumed  equal  to  the  temperature  of  the  material  in  the  bed. 

For  the  shale  bed  previously  defined  with  27%  porosity,  the  void  volume  is 
7  3 

approximately  1.2X10  ft  .  The  total  free  flow  area  is 
A  =  V^H  =  1.66  X  10^  . 

The  heat  required  to  liberate  all  oil  from  the  nuclear  chimney  may  be  determined  by 
q  =  hw  ,  (3) 

where  h  is  the  enthalpy  and  w  the  weight  of  oil  shale,  or 
q  =  1.206  X  10^^  Btu. 


Parameters  of  use  in  determining  overall  heat  transfer  characteristics  are 

listed.  The  effective  heat  transfer  coefficient  h  is  defined  as 

e 


h 


e 


NuX 


(4) 


The  effective  diameter  d  is  assumed  to  be  2  ft. 

e 

The  average  temperature  difference,  AT,  between  the  heat  transporting  fluid  and 
the  bed  solids  is  defined  as 


T^  -  Tg  -  37op 
=  ln[(Tj  -  Tg)/3] 


It  is  assumed  that  no  significant  heat  transfer  occurs  for  temperature  differences  less 
than  37°F. 

The  heating  surface,  S^,  per  foot  of  bed  area  parallel  to  the  direction  of  flow 
within  bed  height,  H,  is  defined  as 

=  [6(1  -  ^)H]/d^  .  (6) 

By  utilizing  these  parameters,  the  heat  transfer  is  defined  as 

Q  =  h  S  AT  .  (7) 

^  e  V 
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Also, 


Q  =  •wpjpC^(Tj  -  Tg)  . 


(8) 


The  average  temperature  difference  AT  does  not  appear  because  its  value  is 
insignificant  compared  to  T^  “  ^s"  obvious  that  the  fluid  velocity  may  be  computed 

from  Eq.  (8).  Based  on  Eqs,  (4)  through  (8),  the  calculations  for  1  through  6  mo  of 
retorting  times  are  made. 


One -Month  (720  hr)  Retort  Time 
The  required  heat  transfer 

Q  =  q/(tA)  (9) 

=  1.206  X  10^^/ (7. 2  X  10^  X  1.66  X  lO'^) 

=  1.01  X  10^  Btu/hr  ft^  . 

Therefore,  the  flow  velocity 
■w  =  Q/[PfC:j(Tj^  " 

=  1.01  X  10^/[54.9X  0.775  (700  -  77)] 

=  3.81  ft/hr. 

The  effective  heat  transfer  coefficient 
hg  =  Q/S^AT 

=  1.01  X  10^/1.585  X  10^  X  1.1  X  10^ 

=  0.58  Btu/ft^  hr  “F  . 

Since 

Nu  =  hg(dg/Xg) 

=  0.58  X  2/6.33  X  lO"^ 

=  18.45, 
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and 


Nu  =  1.25  X  10"^ 


then  the  Reynolds  number 


Re  =  708  . 


By  definition 


Re 


d  G 
e 


Then, 


and 


G  =  Re  |u/d^ 


=  145.5  Ib/ft  hr. 


Wj  =  GA 


=  2,420,000  Ib/hr 


=  330,000  gph. 

12 

To  compute  the  heat  transfer  coefficient,  h^,  between  the  oil  and  the  shale 


h  d 


^  =  2.0  +  0.6  Pr^/^  Re^/^  , 

k  * 

s 


where  the  Prandl  number  is  defined  as: 


Pr  = 


then 


h^  =  1 8. 1  Btu/hr  ft"^ 


(10) 
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If  the  prior  data  are  used  to  compute  temperature  transient  times  in  infinite  slabs, 

13 

the  following  characteristics  are  noted.  According  to  McAdams,  heat  transfer 
problems  on  infinite  slabs  may  be  reduced  to  dimensionless  parauneters  X  and  Y: 

Y  =  (Tj  -  T)/(Tj  -  Tg)  . 

The  symbols  r  and  r^  are  related  according  to  Fig.  3.  If  r^  is  the  half  thickness  of 
the  infinite  oil  shale  slab  and  r  is  the  distance  from  the  midplane,  then  the  parameters 
m  and  n  are  defined  such  that 


and 


n 


The  time  it  would  take  the  temperature  of 
the  slab  to  be  raised  to  690°F  at  its  center 
plane  if  the  slab  is  heated  from  both  sides 
is  now  determined.  If  r  =  0  and  r^  =  1  ft, 
then  n  =  0  and  m  =  0.07.  In  the  computation 


Fig.  3.  Infinite  slab.  Y  =  (700  -  690)/(700  -  77) 

=  0.016. 
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Observing  the  charts  on  pp.  32  to  35  of  McAdams  for  m  =  0.07  and  Y  =  0.016,  it  is 
seen  that 


X  =  1.8  . 


Solving  for  the  time,  t,  it  can  be  seen  that 
t  =  66  hr. 
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When  the  slab  is  heated  from  one  side  such  that 


m 


2  ft 


m  =  0.0358, 


assume 


n  =  0.8  (r  =  1.6  ft) 
Y  =  0.016  ; 


therefore,  it  is  observed  that 
X  =  1.3 

and  the  time  it  takes  the  temperature  to  reach  690®F  at  a  distance  of  1.6  ft  from  the 
heated  edge  is 


t  =  191,5  hr. 


Three-Month  (2160  hr)  Retort  Time 


The  following  are  the  required  and  computed  parameters: 

Q  =  3,37  X  10^  Btu/hr  ft^ 
w^  =  1.27  ft/hr  (oil  velocity) 
hg  =  0.1935  Btu/ft^  hr  “F 
Nu  =  6.11 
Re  =  338.9 
G  =69.6  Ib/ft^  hr 

=  1,148,700  Ib/hr  (156,500  gph) 
h^  =  12.9  Btu/hr  ft^  °F. 

Heating  an  infinite  oil  shale  slab  from  both  sides,  t  =  77.8  hr  to  raise  the  midplane 
to  690®F,  Heating  from  one  side,  t  =  287  hr  to  raise  the  temperature  to  690°F  at  1.6  ft 
from  the  heated  edge. 
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SiX"Month  (4320  hr)  Retort  Time 


The  following  are  the  required  and  computed  parameters; 


Q 

w 

o 

h 

e 

Nu 

Re 

G 

W 

o 


1.68  X  10^  Btu/hr 
0.636  ft/hr  (oil  velocity) 
0.09638  hr  “F 

3.05 
208.9 

43  Ib/ft^  hr 

712,700  Ib/hr  (97,100  gph) 
10.4  Btu/hr  ft^  °F  . 


In  heating  an  infinite  oil  shale  slab  from  both  sides,  the  time  t  =  80.2  hr  to  raise  the 
midplane  to  690°F,  Heating  from  one  side,  t  =  294  hr  to  raise  the  temperature  to 
690°F  at  1.6  ft  from  the  heated  edge. 

In  observing  the  flow  velocity  and  the  heating  transients  of  the  infinite  slab 

problems,  it  is  apparent  that  retorting  times  greater  than  three  months  are  excessive. 

Economics  should  dictate  shorter  retort  times  in  the  range  comparable  to  twice  the 

number  of  hours  to  heat  an  infinite  slab  from  both  sides.  It  may  be  necessary  to  fill 

the  void  spaces  with  oil  at  elevated  temperatures  and  to  maintain  a  flow  predictable  by 
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Darcy  Flow  type  equations. 


Oil  Flow  Patterns 


The  Darcy  Flow  equation  is  of  the  form 


U  =  - 


AC 


+  gp  sm  6 


sec  ’ 


where 


C  = 


5SJ  (1  -  <|>)^ 


(11) 


The  angle  6  is  taken  with  respect  to  the  horizontal  direction.  For  flow  downward, 

0  =  90°.  It  can  then  be  seen  that  the  sin  0  =  -1,  If  no  flow  takes  place  with  respect  to 
the  referenced  vertical,  the  fluid  is  in  static  equilibrium  and  obviously 
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^  -  gp  =  0  . 

K  the  bottom  of  a  porous  medium  is  vented  and  the  rate  of  flow  is  influenced  only 
by  gravity,  then. 


AP 

H 


gP  <  0  , 


due  to  a  small  AP.  The  maximum  rate  of  flow  due  to  gravity  takes  place  when  AP  -►  0. 
By  inducing  a  slight  vacuum  at  the  bottom  of  a  porous  medium  with  respect  to  the 
atmosphere,  AP  can  be  influenced  to  the  extent  that  the  flow  rate  will  increase.  The 
Darcy  Flow  equation  is  applicable  with  systems  of  variable  porosity;  i.  e. ,  elements  of 
height  may  be  considered  if  the  porosity  within  a  defined  element  is  assumed  constant. 


Temperature  Profiles  in  Shale  Bed 


The  solution  to  the  problem  of  the  temperature  profile  in  a  porous  medium  due  to 
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a  hot  flowing  liquid  was  approached  from  the  formulations  established  by  Schumann, 

A  right  circular  cylinder  consisting  of  crushed  material  at  some  uniform  temperature 
was  assumed.  A  fluid  at  some  higher  uniform  temperature  passes  through  the  crashed 
material  at  a  uniform  rate  of  flow.  The  sides  of  the  cylinder  are  assumed  impervious 
to  the  fluid.  The  distribution  of  temperature  in  the  crushed  material  and  the  fluid  for 
all  times  is  found  based  on  the  assumption  that 

1.  The  lumps  of  material  are  so  small  or  have  such  high  thermal  diffusivity  that 
any  given  lump  may  be  considered  at  a  uniform  temperature  at  any  instant. 

2.  Compared  to  the  transfer  of  heat  from  fluid  to  solid,  the  transfer  of  heat  by 
conduction  in  the  fluid  itself  or  in  the  solid  itself  is  small  and  may  be  neglected. 

3.  The  rate  of  heat  transfer  from  fluid  to  solid  at  any  point  is  proportional  to  the 
average  difference  in  temperature  between  fluid  and  solid  at  that  point. 

4.  Change  in  volume  of  fluid  and  solid  due  to  change  in  temperature  may  be 
neglected. 

5.  The  thermal  constants  are  independent  of  the  temperature. 

In  the  application  of  Schumann’s  efforts  to  the  oil  shale  problem,  it  became 

apparent  that  the  curves  generated  by  him  did  not  extend  far  enough  to  give  solutions 

for  the  ratio  of  the  fluid  to  the  solid  temperatures.  It  was  necessary  to  investigate  the 
1 6 

work  of  Furnas.  These  curves  also  were  inadequate;  therefore,  it  was  necessary  to 
derive  new  equations  the  solution  of  which  involved  Bessel  Functions  with  large 
arguments. 
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According  to  Schumann,  the  equation  expressing  the  temperature  ratio  between  a 
solid  and  a  fluid  may  be  written  as 


^  =  1  -  e-Y-Z 


^  y"m^(YZ)  , 


(12) 


n=0 


where  the  arguments  Y  and  Z  are  defined  as 


Y  =  h^x/(CjVj<t>) 


and 


Cs  >1  -  « 


The  symbol  x  is  the  distance  traveled  by  the  fluid  front  into  the  porous  medium. 
If  the  following  transformations  are  made: 

Y  +  Z  =  $ 


YZ  =  § 


then  Eq,  (12)  may  be  rewritten  as 


=  1  -  e-® 


(13) 


n=0 


17 

By  definition 


8  M.(l) 
M  (I)  = - 2 - 

n'  ggn 


(14) 


where 


2  3 

mJ?)  =  1  +  c  +  -L^  +-5 

^  {2'.r  {3\r 
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In  expanding  Eq,  (13), 


-  =  1  -  e'®  [mq(5)  +  YMj(?)  +  Y^M2(?)  +  . . .  j  , 


but  utilizing  Eq.  (14) 


T,  .*  [  3m.(|)  ,  a^M. 


The  function  M^CS)  is  usually  defined 


Mq(§)  =  Iq  (2^/f )  s  y  a)  , 


where 


u  -  2n/? 

and  Iq((t)  is  a  modified  Bessel  Function,  Keeping  in  mind  that 

-rV  -  V»> . 


Eq.  (16)  may  be  rewritten 


ya) 


+  Y 


w 


9a 


+  Y 


[inf 


8»£2!V!^U  1 

8?  352  8c  I  • 


Also, 


8c7  ^  2 
■  <7* 


(15) 


(16) 


(17) 
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If  the  modified  Bessel  Functions  are  expanded  into  infinite  series  and  all  higher 
order  terms  neglected  for  cr  »  1,  then 


T 


=  1  -  — - 

■f--(-) 

n/  27ra 

CJ  \\c 

7/  cr  \g/ 

)  J 

or 


T 

_s«  1 

^o  ^[2^0 


1+2X  + 

CT 


<! 

'  1 

1 

N  -  JL  ’ 

1  ‘’1 

>+  .  .  . 

=  1 


“^a 


n/  27rcr 


1  + 


(18) 


/2Yy 

<  1  such  that  Y  <  Z,  then  this  solution  of  the  temperature  ratio  is  good 

only  for  time 


t>^ 


1  + 


s  (1  "  <^) 

Cj) 


(19) 


therefore. 


T„  ^-4>+or 

=  1 


(20) 


and 


c|)  =  Y  +  Z 

CT  =  2^^YZ  . 


Design  of  Heat  Exchangers 


Heat  exchanger  design  techniques  are  similar  to  those  presented  for  Sodium  (Na) 


18 


and  sodium-potassium  (Na-K)  liquid  metal  systems  by  Fraas  and  Ozisik,  In  the 


-16- 


proposed  triple  cycle  heat  exchanger,  oil  is  the  third  cycle  coolant,  but  the  underground 
oil  shale  represents  the  heat  sink.  In  each  heat  exchanger  the  fluid  flows  in  a  parallel 
manner. 

On  the  initial  pass  it  is  assumed  that  the  oil  enters  the  heat  exchanger  at  a 
temperature  equal  to  the  initial  oil  shale  temperature  (77°F).  The  oil  will  be  raised 
to  a  temperature  of  700'^F  before  being  injected  into  the  cool  oil  shale  chimney.  From 
the  equation 

Q  .  W„C„AT, 

the  heat  flow  to  the  oil  for  a  30 -day  retorting  system  is 
Q  =  1,17  X  10®  Btu/hr. 

The  parameters  governing  the  Na-K  secondary  heat  exchange  cycle  are  listed  below. 


Sodium-Potassium  Secondary  Heat  Exchange  Cycle 
The  following  are  the  required  and  computed  parameters: 


Tl 

=  795°F 

^2 

=  945°P 

^NK 

=  0.2503  Btu/lb°F 

%K 

=  16  Btu/lb  ft  “F 

H 

=  0.494  Ib/hr  ft 

=  51  Ib/ft^ 

Pr 

=  0.00842 

PrO.4 

=  0.14775 

i,  d. 

=  0.305  in. 

o.  d. 

=  0.375  in. 

The  subscript  NK  means  sodium-potassium  (Na-K). 

For  the  oil  heating  cycle,  the  Na-K  flows  inside  small  tubes,  while  the  oil  flows 
against  the  outer  surfaces  of  the  tubes.  An  expression  for  the  required  number  of  tubes 
is 


W 


n  = 


NK 


6.31 


fL 


-]l/2 


(21) 
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If  the  heat  load  between  the  oil  and  Na-K  is  the  same,  then  for  the  oil: 


^o  "  ^NK 


Ww  =  r - -  =  8,61  X  10^  Ib/sec. 

NK  Cj^k^-^NK 


Assume  f  =  0,03,  L  =  10  ft,  =  51  Ib/ft^  and  AP  =  1440  Ib/ft^  then 


n  =  32,200  tubes. 

The  flow  rate  per  tube  is  then 

W  =  W  /n 
NK  NK' 

=  0,2674  Ib/sec, 

The  tube  inside  diameter  (i,  d, )  is  assumed  to  be  (D  =  0,305)  in, ;  therefore,  its 
cross  sectional  area 

A  =  0,0731  in^. 

The  mass  flow  rate  per  tube  is 


W 


NK  A 


=  526  Ib/sec  ft^. 


and  the  Reynolds  number 


Re  =  =  97,500  , 


The  heat  transfer  coefficient  for  Na-K  is  computed  from  the  expression 


hj^j^  =  0,02  • 


(22) 


From  the  known  parameters 


h™^  =  55,200  Btu/hr  ft'^  °F  . 
JNJx 
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9 

On  the  side  of  the  oil,  Q  =  1.17  X  10  Btu/hr.  The  total  oil  flow  rate  through 
the  exchanger  is 

W  =  673  Ib/sec. 
o  ' 

Outside  the  tube,  the  flow  rate  per  tube 
=  0.0209  Ib/sec, 

I  2 

and  the  mass  flow  rate  =  37.8  Ib/sec  ft  at  a  Reynolds  number 
Re  =  10,370. 

The  Prandl  number  computes  to  be  (Pr  =  5.04);  therefore,  the  heat  transfer 
coefficient 

=  177  Btu/hr  ft^  “F, 

A  columbium  tube  wall  is  assumed  between  the  two  media: 

=  34  Btu/hr  ft  °F  . 

The  wall  conductance  (heat  transfer  coefficient) 
h  =  k  /thk,  thk  =  0.035  in. 

=  11,650  Btu/hr  ft^  °F. 

The  overall  heat  transfer  coefficient  h  is  derived: 

u 

h  h,^„  h  h  (i.  d./o.  d.) 

u  NK  w  o  ^  ' 

The  tube  outside  diameter  is  0.375  in.  and  the  inside  diameter  is  0.305  in. 
h^  =  138.5  Btu/hr  ft^  “F  . 

The  total  tube  external  surface  required  is  calculated  by  equating  the  power  from 
the  Na-K  to  the  oil: 

Q  =  Ah  AT 
^  u 
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where 


AT 


NK 


In  [(T 


-  T  -  37 

-°V/3] 


NK 


Remembering  that  -  945°F,  -  77®F;  then 

AT  =  147°F. 


(24) 


Since 

Q  =  1.17  X  10®  Btu/hr, 

A  =  5,74  X  10^  ft^. 

The  external  surface  per  tube  is 

a  =  A/n  =  1.78  ft^/tube. 

The  tube  length 

L  =  A/C, 

where  C  is  the  circumference  of  the  outside  diameter,  or 
L  =  18.15  ft. 

For  the  Na-K  pressure  drop 

=%K^NK^/<2gpD).  (25) 

3  2 

Assume  =  0.03;  therefore,  AP  =  1,885  X  10  Ib/ft  (13  psi)  which  is  fairly  close; 
JNK  , 

therefore,  3,22  X  10  tubes  are  a  good  approximation.  For  the  oil,  if  f^  -  0.035, 

AP  =  8,23  Ib/ft^. 
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Sodium  Primary  Heat  Exchange  Cycle 


The  following  are  the  required  and  computed  parameters: 

Tl  =  950°F 

T2  =  1050°F 

=  0.302  Btu/lb  °F 

=  37,2  Btu/lb  ft  °F 

ju  =0.56  Ib/hr  ft 

=51,2  Ib/ft^ 

Pr  =  0.00455 
0  4 

Pr  '  =  0.1155 

AP  =  1440  Ib/ft^  (allowed) 
i,  d.  =  0.305  in. 

Q  =  1.17  X  10^  Btu/hr 

The  subscript  N  refers  to  sodium  (Na). 

Employing  the  same  expressions  used  for  the  secondary  cycle,  the  following 
are  the  parametric  requirements  computed: 

=  10,750  Ib/sec 
n  =  33,600  tubes 

=  0.3212  Ib/sec  tube  (inside) 

Wnk  "  0.257  Ib/sec  tube  (outside) 

The  mass  flow  rates  are: 

2 

=  633  Ib/sec  ft  (inside) 

Gnk  "  Ib/sec  ft^  (outside) 

Re^  =  99,000 
Re,,„=  1.07  X  10^ 

X\ 

The  heat  transfer  coefficient  values  are: 

=  33,700  Btu/hr  ft^  “F 
h^K  =  30,200  Btu/hr  ft^  °F. 

Between  the  Na  and  the  Na-K,  the  average 

at  =  48®F. 
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Assume  a  columbium  tube  wall  at  lOOO'^F,  Then 


k^/p  =  0.067 

-  35.8  Btu/hr  ft  °F 
thk  =  0.035  in. 

h  =  12,300  Btu/hr  ft^  °F. 

w 


The  overall  heat  transfer  coefficient 
=  2420  Btu/hr  ft^  °F. 

The  total  required  external  surface 
A  =  10,000  ft^. 

The  tube  length 

L  =  4.82  ft. 

The  pressure  drop  on  both  sides  of  the  tube  becomes 

APpg  =7  61  Ib/ft^  (inside) 

=  824  Ib/ft^  (outside). 

Therefore,  the  number  of  tubes  is  reasonable  since  the  pressure  drops  are  less  than 
those  assumed. 

For  all  designs,  the  tube  spacing  should  be  square  and  at  least  1.25  times  the 
tube  outside  diameter  or  larger.  Heat  exchanger  designs  should  be  comparable  to  those 
developed  by  the  AEC  for  liquid  metals  and  molten  salts  for  best  performance. 


Nuclear  Reactor  Heat  Transfer 


In  the  computation  of  fuel  element  temperatures,  it  is  assumed  that  the  lengthwise 
heat  distribution  is  generated  in  a  cosine  pattern.  The  maximum  fuel  element  temperature 
in  Eq.  (26)  depends  upon  three  basic  reactor  characteristics;  i.  e. ,  core  coolant  inlet 

temperature  T..,  coolant  temperature  rise  AT,  and  the  average  coolant  film  temperature 

^  n  19 
drop  0^. 


-22- 


The  expression  for  the  maximum  fuel  element  surface  temperature 


AT 

T  =  T.  +  ^  + 

m  1  AT  2 


1/2  J 


X  rr.2  .  -^2  2  ^2 

+Fe' 


(26) 


where  the  average  temperature  rise  in  the  core 


AT  =  0.008 


S  L 

L^P  ^  ^  » 


0.357 


n 


0.0714  Q 

jU _  ^  _  _ 

_0.429  C  qO.643  _0.929 

Dg  p  /3  D 


(27) 


and  the  average  coolant  film  temperature  drop  in  the  core 


CQ  o  r.0-2  ArpO.8  1.2  ^  0.4  0.4 

58.2  AT  De  Cp  p 

t^1.4  ,0.6  qO.2  078 

D  (Lq/D)  /3  n 


(28) 


The  hot  channel  factor  accounting  for  nonuniform  racdial  heat  generation  and 
average  dimensional  deviations  is  chosen  to  be  ~  1,5,  For  a  coolant  film 
temperature  drop  small  compared  to  the  average  temperature  rise  in  the  core  when 

AT  =  lOO^F  , 

the  maximum  fuel  element  surface  temperature  is  approximately 

T  «  1200°F. 
m 


Sodium  Graphite  Nuclear  Reactor 

The  following  are  the  required  parameters: 

T^  =  950‘^F 
T^  =  1050°F 
T  =  1200°F 

w  n 

p  =  51.2  Ib/ft 

The  subscript  w  refers  to  the  reactor  coolant  flow  walls. 

The  sodium  flow  rate  through  the  reactor  core  is  the  same  as  through  the  primary 
heat  exchanger: 


Assuming  a  coolant  passage  inside  diameter  of  1  in. ,  an  allowable  pressure  drop  of 
10  psi,  and  a  10 -ft  core  length,  the  following  are  the  resulting  parameters  computed: 

n  =  5440  tubes 
=  1.98  Ib/sec 
G^  =  363  Ib/sec  ft^ 

Re^^  =  195,000 

=  1,74  X  10“^  Btu/hr  °F 

Assuming  a  wall  temperature  of  1200°F  and  the  coolant  inlet  temperature  of  950®F,  the 
average  temperature  difference  between  the  walls  and  the  bulk  fluid  becomes 

AT  =  57°F 

h^  =  1520  Btu/hr  ft^  °F. 

4  2 

The  coolant  total  free  flow  area  A  =  1.35  X  10  ft  .  Therefore, 
a  =  2.48  ft^/ tube. 


The  computed  L  =  9.5  ft,  and  the  AP  for  9.5  ft  is  144  psf,  but  the  AP  for  15  ft  is  216  psf. 
For  a  square  core  where  the  diameter  and  the  length  equal  9.5  ft,  the  volume  fraction 
of  sodium  in  the  core  causes  a  high  fuel  inventory.  It  is  necessary  to  assume  a  core 
diameter  and  length  of  15  ft,  each.  For  a  1-in.  diam  flow  passage  the  total  sodium 
frontal  area  is  computed  to  be 

=  29.648  ft^. 

The  reactor  frontal  area 

A  =  176.7  ft^. 
r 

The  reactor  core  volume 

V  =  2650.5  ft^  =  7.506  X  10^  cm^, 
r 

The  sodium  void  fraction  equals  16.77%. 
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For  a  100°F  temperature  rise  in  the  core 


10,750  X  3600  X  0.302  X  100 
5440  X  3413 

=  63  kw/channel 

Therefore,  the  total  reactor  power  is 

Q  =  340  MW. 
r 


Reactor  Criticality  Calculations 


The  volume  fraction  of  the  core  materials  are  as  indicated  in  Table  IL 


Table  II.  Reactor  composition 


%  Vol 

Element 

Volume 

(cm^) 

Weight^ 

(g) 

Core  density 
(g/cm^) 

barns 

3 

Atoms/cm 

""a 

0.1677 

Na 

1,26  X  10^ 

1.17  X  10^ 

0.1555 

0.4 

4 

4.07  X  10^^ 

0.0310 

ss 

2,33  X  lo"^ 

1.86  X  10^ 

0.2480 

2.4 

11 

2.67  X  10^^ 

0.0500 

Nat.  U 

3.75  X  lo"^ 

7.13  X  10^ 

0.9500 

2,40 X  10^^ 

— 

U235 

— 

— 

—  687 

9 

— 

— 

U238 

— 

— 

— 

2.75 

8.3 

— 

0.7513 

Graphite 

5.64  X  10^ 

o 

o 

X 

o 

cc 

1.3345 

0.0045 

4.8 

6.70 X  10^^ 

1,0000 

7.51  X  10^ 

2.016  X  10® 

^Total  weight  without  shield  (4.445  X  10^  Ib). 


For  natural  uranium  1,111%  enriched  with  U235,  the  ratio  of  the  atomic  densities 
^238^^235  solving  the  expression 

^  '^a235  ^^238^  •'^235^  °'a238 

au  1  + 

the  thermal  absorption  cross  section  for  the  enriched  natural  uranium  is  computed  to 
be  =  10.27  barns.  The  macroscopic  absorption  cross  section  =  0,025  cm  ^  for 
the  enriched  uranium  where  S  =  Na.  The  scattering  cross  section  is  computed  with  a  . 
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similar  expression  such  that  the  macroscopic  scattering  cross  section  =  0.0199  cm 

In  summing  the  macroscopic  cross  sections  for  the  absorption  and  the  scattering 

process  separately,  the  total  absorption  and  scattering  macroscopic  cross  sections  are 

2  =  0.0339  and  S  =  0.39716,  respectively, 
a,  s 

From  the  parameters  known,  a  thermal  utilization  factor  f  may  be  computed, 

where 


f  =  2  /  2^ 

au'  a 

=  0.725 

represents  the  percent  of  neutrons  utilized  to  carry  on  fissions. 

The  number  of  fast  fission  neutrons  produced  per  thermal  neutron  absorbed  in 

.  20 

uranium  (rj)  is  computed  from  the  expression 

_ _ 2.5 

^  1.184  +  0.00474  (N238/N235)  ■ 

For  this  case,  n  =  1.552.  The  number  of  second  generation  fissions  in  U235  per  fission 
of  U235  from  the  first  generation  in  a  reactor  of  infinite  extent  is  computed  as  the 
multiplication  factor  defined  by 

rif 

K  =  1.125  . 

00 

The  diffusion  coefficient  D  may  be  computed  from  the  expression 

°  "  323  (1^  -  Pq)  ’ 

where  =  2/3A.  Since  this  is  a  homogenous  reactor,  A  represents  the  weighted 
average  of  the  atomic  weights  in  the  core: 

A  =  21.5091 
=  0.030994 
D  =  0.86612  cm. 
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The  square  of  the  diffusion  length  L  is  computed  from 


=  D/S 
'  a 

=  25.4988  cm^  . 


For  a  pure  graphite  moderator,  the  slowing-down  length  L  =18.7  cm.  If  only 

s 

the  moderator  influences  a  neutron  lifetime  in  the  core  materials,  then  the  modified 
slowing-down  length 


density  of  the  graphite 


L  =18  7 _  _ 

s  ‘  density  of  graphite  per  unit  volume  of  core 


=  22.406  cm. 


By  definition. 


T  '2 

Lg  =T  , 


,2.  ... 


the  Fermi  age. 

For  a  reactor  of  the  given  dimensions,  the  buckling  term  (B"^)  is  computed  by 
the  expression 


2  -4 

B  =  1.5456  X  10 


The  effective  multiplication  factor  due  to  a  finite  geometry  cam  now  be  computed: 


xr  =  ^ _ 

I  +  lV 

=  1.0369 

The  excess  multiplication,  K  =  K  « «  -  1  =  0.0369,  will  determine  the  reactor  lifetime. 

IQ  3  fi 

For  a  reactor  where  =  1,  =  2.45738807  X  10  atoms/ cm  (7.2  X  10  g),  the 

thermal  neutron  flux 


1.2  X  10^^  P 

1 

Wa^ 


The  excess  reactivity  which  affects  the  reactor  lifetime  was  an  arbitrarily  built-in 
factor. 
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•where  W  is  the  weight  of  U235  in  grams  and  the  reactor  power  in  watts.  For 

12  2 

ip  =  1.03  X  10  neutrons/cm  -sec 


and 

Sj  =  0.0135  cm"^  , 

the  average  number  of  fissions  per  second  ^is  determined  by 

19  f 

For  =  1,  the  fission  rate  ^  =  1.04  X  10  fissions/sec.  The  lifetime  of  this 
reactor  with  an  arbitrarily  selected  3.69%  reactivity  is  calculated  to  be 

t  =  N235  V/^sec 

t  =  1.8947  X  10^  sec. 

Since  there  are  3.145  X  10  sec/year,  therefore,  t  =  63  years. 

The  reactor  lifetime  may  be  reduced  to  any  desired  number  of  years  by  reducing 
the  fuel  inventory.  Also,  for  the  sake  of  mobility,  the  reactor  size  may  be  reduced, 
governed  by  the  allowed  coolant  pressure  drop  across  the  core.  Any  reductions  cause 
favorable  sa-vings  in  the  overall  system  weight. 


Theoretical  Pumping  Power 


Computations  on  pumping  power  make  no  allowance  for  pump  efficiency.  The 
following  equation  expresses  pumping  power: 

Pp  =  (AP)  (A^)  (v^)/(33,000)  (60) 

For  the  primary  cycle,  the  following  parameters 

Arj,  =  17.05  ft^ 

AP  =761  Ib/ft^ 

=  4,44  X  10^  ft/hr 

yield  a  pumping  power  of 

P  =  290.5  hp  (minimum). 

P 
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For  the  secondary  cycle,  the  following  parameters 


=  16.35  ft^ 

AP  =  1885  Ib/ft^ 

=  3.72  X  10^  ft/hr 

yield  a  pumping  power  of 

P  =  578  hp. 

P 

For  the  oil  cycle,  the  following  parameters 

=  2.4228  X  10®  Ib/hr 

H  =  3750  ft  (pumping  height) 

P 

yield  a  pumping  power  of 

Pp  =  H  /(33,000)  (60) 

=  4590  hp. 

The  pumping  power  requirement  through  the  heat  exchanger  is  negligible  compared  to 
the  pumping  requirement  from  the  bottom  of  the  chimney. 


Factors  Affecting  System  Economics 


An  economical  appraisal  of  this  system  is  beyond  the  scope  of  this  report.  Some 
21  22 

of  the  factors  ^  affecting  the  economics  of  this  system  are  commented  upon.  Con¬ 
sideration  must  be  given  to  the  reactor  system,  heat  exchanger,  pipe  lines,  land 
leasing,  roads,  interest  and  amortization  rates,  operational  cost  over  specific  time 
spans,  moving  of  the  reactor,  nuclear  explosive  services  and  manpower. 

The  fact  that  66  hr  is  required  to  retort  a  2-ft  thick  oil  shale  slab  at  a  background 
temperature  of  700®F  suggests  that  a  minimum  of  10  nuclear  chimneys  may  be  retorted 

in  4  wk.  This  is  equivalent  to  132  chimneys  per  year.  For  an  oil  shale  field  of  a 

0 

28-gpt  yield,  the  oil  shale  tonnage  for  the  given  chimney  is  2,5  X  10  tons.  The  oil 

7  6 

yield  would  be  7  X  10  gal  or  2,22  X  10  barrels  of  oil  per  chimney. 

The  lifetime  of  a  nuclear  reactor  will  affect  the  economics.  A  nuclear  reactor 
can  be  designed  for  any  lifetime.  Coupled  to  the  lifetime  is  the  reactor  power 


-29- 


requirement  of  diminishing  output  over  the  retorting  time  of  a  single  chimney.  As  the 
oil  shale  approaches  retorting  temperature,  the  reactor  power  requirement  falls  off. 

Another  factor  affecting  the  economics  is  the  residual  heat  in  a  nuclear  chimney. 
The  residual  heat  in  the  oil  shale  is  a  function  of  time  after  chimney  formation.  For 
example,  the  maximum  observed  temperature  11  mo  after  detonation  in  grandiorite 
was  88°C  for  the  Hardhat  experiment. 

The  nuclear  reactor  power  requirement  will  be  governed  by  the  equilibrium 
temperature  in  a  nuclear  chimney.  Nevertheless,  the  size  (weight)  of  a  nuclear 
reactor  and  its  mobile  unit  will  be  governed  by  the  required  energy  the  reactor  must 
supply  to  take  the  oil  shale  to  retorting  temperatures. 


System  Advantages  and  Disadvantages 


System  Advantages 

This  system  has  the  following  advantages: 

1.  The  maximum  temperature  of  retorting  is  controllable. 

2.  The  retorting  time  is  controllable. 

3.  When  a  hydrocarbon  is  employed,  no  chemical  reaction  is  expected  between 
the  hydrocarbon  and  the  oil  shale. 

4.  Residue  from  the  oil  shale  retorting  process  is  left  in  a  favorable  condition 
for  additional  reclamation  processes.  No  new  compounds  have  been  formed  as  a 
result  of  the  retorting  process. 

5.  The  nuclear  chimney  walls  give  up  oil  to  the  retorting  process.  Because  of 
this,  oil  yield  can  exceed  the  calculated  yield  of  the  broken  oil  shale  within  the  chimney, 

6.  Oil  losses  are  kept  at  a  minimum  by  constant  temperature  medium  retorting. 

7.  Oil  retorted  from  residual  heat  in  neighboring  chimneys  may  be  used  to 
start  the  retorting  process. 

8.  Only  one  externally  supplied  oil  reserve  is  required  to  start  the  retorting 
process  over  the  life  of  the  nuclear  reactor, 

9.  By  centrally  locating  a  mobile  retorting  system  among  several  nuclear 
chimneys,  the  cycling  of  oil  through  the  heat  exchanger  need  not  terminate  with  the 
emptying  of  a  single  nuclear  chimney. 

10.  The  residual  heat  in  a  nuclear  chimney  remaining  from  the  retorting  process 
may  be  used  to  generate  electrical  power  by  steam  generators  or  MHD  power  con¬ 
version  devices. 

11.  The  rate  of  oil  extraction  is  controllable. 
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12.  The  minimum  time  for  oil  shale  to  be  brought  up  to  700°F  is  approximately 
66  hr  (2.75  days).  Two  chimneys  should  be  retorted  completely  in  1  wk. 

13.  A  diesel-driven  mobile  retorting  unit  with  appropriate  electrical  generators 
could  supply  the  necessary  power  for  all  pumps. 

14.  Fluids  other  than  hydrocarbons  may  be  used  as  heat  transporting  media. 

15.  An  alternate  to  nuclear  power  is  the  direct  application  of  heater  coils  for 
heating  the  heat  transporting  medium. 


System  Disadvantages 

The  following  are  disadvantages  of  the  system: 

1.  The  mobile  oil  shale  retorting  unit  is  heavy. 

2.  Building  more  than  one  mobile  oil  shale  retorting  device  would  require  an 
outlay  by  the  Federal  Government  or  by  a  combination  of  several  industrial  firms. 
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Appendix  A 
List  of  Symbols 


D  ,  d 
e  e 


area  (ft^) 

weighted  average  of  atomic  weights  in  core 

total  sodium  frontal  area 

reaction  frontal  area 

total  flow  area  (ft^) 

_2 

buckling  term  (cm  ) 

circumference  of  outside  diameter 

constant  to  be  empirically  derived 

specific  heat  to  fluid  at  constant  pressure  (Btu/lb  °F) 

specific  heat  of  sodium  at  constant  pressure 

specific  heat  of  sodium-potassium  at  constant  pressure 

specific  heat  of  oil  at  constant  pressure 

specific  heat  at  constant  pressure 

specific  heat  of  solid  at  constant  pressure 

diffusion  coefficient 

reactor  core  diameter  (ft) 

equivalent  (effective)  diameter  (ft) 

hot  channel  factor  accounting  for  nonuniform  radial  heat  generation  and 
average  dimensional  deviations 

thermal  utilization  factor 

friction  factor 

sodium  friction  factor 

sodium-potassium  friction  factor 

hot  channel  factor  applied  to  coolant  film  temperature  drop  for  nonuniform 
heat  generation  with  average  and  local  off -design  deviations 

ratio  of  system  pressure  drop  to  core  and  reactor  vessel  pressure  drop 

ratio  of  core  and  reactor  vessel  pressure  drop  to  core  friction  pressure 
drop 

mass  flow  (Ib/hr  ft^) 

mass  flow  for  sodium 

mass  flow  for  sodium-potassium 

acceleration  of  gravity  (ft/sec  ) 

pumping  height  (ft) 

bed  height 

height  of  chimney  (m) 
enthalpy 
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h  depth  of  burial  of  nuclear  device  (m) 

heat  transfer  coefficient  between  oil  and  shale 
h^  effective  heat  transfer  coefficient 
heat  transfer  coefficient  of  sodium 
^NK  transfer  coefficient  of  sodium-potassium 

h^  heat  transfer  coefficient  of  oil 

h^  overall  heat  transfer  coefficient 

♦  h  wall  heat  transfer  coefficient 

w 

Iq{(j)  modified  Bessel  function 

K  empirically  derived  constant 

thermal  conductivity  of  a  fluid  (Btu/ft  hr  °F) 

thermal  conductivity  of  sodium  (Btu/ft  hr  ®F) 

^NK  thermal  conductivity  of  sodium-potassium  (Btu/ft  hr  °F) 

k^  wall  thermal  conductivity  (Btu/ft  hr  ®F) 

kg  thermal  conductivity  of  a  solid  (Btu/ft  hr  °F) 

Keff  effective  multiplication  factor 

excess  multiplication  factor 
ex 

multiplication  factor 
L  length  (ft) 

L  slowing  down  length  (cm) 

s 

In  natural  logarithm 

N  atomic  density 

Nu  Nusselt  number 

m  parameter 

*  n  number 

n  parameter 

P  pressure  (psf) 

P  Prandl  number 

r 

P^  reactor  power  (W) 

Pj^  pressure  for  sodium  (psf) 

Pnk  pressure  for  sodium-potassium  (psf) 

P  pumping  power  (hp) 

P 

Q  heat  transfer 

heat  transfer  for  sodium 
Qnk  heat  transfer  for  sodium-potassium 

thermal  output  from  reactor  (kW) 
q  heat  required  to  liberate  all  oil  from  nuclear  chimney  (Btu) 

R  cavity  radius  (m) 

Re  Reynolds  number 

^  Reynolds  number  for  sodium 

ROnk  Reynolds  number  for  sodium-potassium 
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r  distance  from  the  midplane  of  an  infinite  slab  (ft) 

r  half  thickness  of  an  infinite  slab  (ft) 

m  2 

heating  surface  (ft  ) 

T  temperature 

maximum  fuel  element  surface  temperature  (°F) 

T  temperature  of  reactor  coolant  flow  walls  (°F) 

w 

inlet  temperature  (°F) 

T2  outlet  temperature  (®F) 

fluid  temperature  (°F) 

initial  fluid  temperature  (°F) 

T  temperature  of  a  solid  (/F) 

s 

t  time 

V  volumetric  flow  rate  (ft  /  sec) 

3 

V  volume  of  cylinder  (ft  ) 

3 

Vy  void  volume  (ft  ) 

flow  velocity 

V  reactor  core  volume 

r 

W  weight  of  (g) 

W  yield  (kt) 

fluid  mass  flow  rate  (lb/ hr) 

mass  flow  rate  of  sodium  (lb/ hr) 

mass  flow  rate  of  sodium -potassium  (lb/ hr) 

mass  flow  rate  of  oil  (lb/ hr) 

w  weight  of  oil  shale  (lb) 

w^  flow  velocity 

w^  oil  flow  velocity 

/3  coolant  volume  fraction 

2 

AP  pressure  change  (Ib/ft  ) 

pressure  change  for  sodium 

APnk  pressure  change  for  sodium-potassium 

AT  average  temperature  difference  (°F) 

AT^  average  temperature  difference  for  sodium 

ATt^tt^  average  temperature  difference  for  sodium-potassium 

NK. 

T}  number  of  fast  fission  neutrons  produced  per  thermal  neutron  absorbed 
in  uranium 

6  average  coolant  film  temperature  change  (°F) 

a 

X  molecular  specific  heat  (Btu/lb  ®F) 

/u  viscosity  (lb/ hr  ft) 

fj,^  fluid  viscosity 

p  average  of  overburden  density  (g/  cm^) 
p.  fluid  density  (Ib/ft^) 
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9 


Pn 

^NK 


CTf 


T 

<!> 

X 


3 

sodium  density  (Ib/ft  ) 

3 

sodium-potassium  density  (lb/ ft  ) 

3 

coolant  density  at  the  pump  (lb/ ft  ) 

mean  coolant  density  in  the  reactor  core  (lb/ ft  ) 

density  of  a  solid 

macroscopic  fission  cross  section  (cm~^) 
macroscopic  absorption  cross  section  (cm”^) 
macroscopic  scattering  cross  section  (cm~^) 
total  macroscopic  absorption  cross  section  (cm"^) 
total  macroscopic  scattering  cross  section  (cm*^) 
microscopic  absorption  cross  section  (barns) 

thermal  absorption  cross  section  for  enriched  natural  uranium  (barns) 

microscopic  fission  cross  section  (cm  ^) 

microscopic  scattering  cross  section  (barns) 

2 

the  Fermi  age  (cm  ) 

%  porosity 

distance  traveled  by  fluid  front  into  porous  medium  (ft) 
average  number  of  fissions  per  second 
thermal  neutron  flux 


f 
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